TITLE 
TYPE 

ACCURACY 
DURATION 



METHOD OF USE 



PUNCHING OF THE DATA 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE M 18 - 213 



One-Step Automatic Eigenvalue -Eigenvector Program 

Entire Program 

Usually 10 or 11 decimal places 

(a) 20 seconds to input program 

2 2 

(b) n /250 to n /kO seconds (depending on number of 

digits ) to input elements of matrix 

(c) n /200 seconds per iteration. Most matrices re- 
quire from k iterations (n = 3) to 7 iterations (n = ho) 
for convergence. 

(d) n /2 seconds to punch results. 

(1) Read program tape into memory in usual way. If 
no read in error has "been committed machine will stop 
on 20 (03L) l6 . 

(2) Place data tape in reader and restart machine. 

(3) After computation has been completed, results will 
be punched out for printing. If only the eigenvalues 
are computed, then they are printed in a single column. 
If the eigenvalues and eigenvectors are computed, then 
following each eigenvalue will be the corresponding 
eigenvector . 

(k) The machine stops on a 20 (03L) 1 ^. A new problem 
can be begun by repeating step 2 . 
To compute the eigenvalues or the eigenvalues and 
eigenvectors of the real symmetric matrix 



A. . 



1 = 1, . . .n; 



j = 1, . . . n proceed as follows : 



(a ) Scale the matrix so that 



n 



i,j=l 



A . < 1/2 . 



(b) Punch the elements of (A. . ) for i > j, row by 
row, as a sign digit followed by up to twelve decimal 
digits . 



CAPACITY 



INTERNAL CHECKS 





LOCATION 


Decim? 


al 


Sexadecimal 


908 


38N 


8 




008 


23 




017 


33 




021 


71 




0^7 


79 




O^L 


168 




0K8 



- 2 - 

('c ) Element A is followed by an N or J . If followed 

v ' nn tf 

by an N, then the eigenvalues and eigenvectors are com- 
puted . If followed by a J, then only the eigenvalues 
are computed . 

(d) The last character is followed by a sexadecimal 
character p which determines the number of decimal 
digits to be printed- The character p can assume the 
values 1,2,3,..., 9, K, or S where K = 10 and S = 11. 

All the eigenvalues and eigenvectors can be found for a 
matrix of order 23° However, all the eigenvalues of a 
matrix of order k-0 can be computed. 

During the input and operation of the program a number 
of checks are made. If they fail, the machine stops. 
A list of the locations and reasons for failure is 
given below 

REASON FOR FAILURE 



Sum check on program tape fails 

Data tape fails to have elements 

Drum failure 

Drum failure 

Drum failure 

Drum failure 

Arithmetic failure in Routine M-4. 

This can be caused by incorrect scaling. 

This program is essentially a combination of M-^-, P-2, 
N-2, and X-80 It replaces .M-fe. 



DAT E Re; April 3, 1957 
PROGRAMMED BY Gene Golub 
APPROVED BY 



ti l (iene LrOiup 



LOCATION 


ORDER 
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Library 


r Routine X-1 

00 3K 
00 F 
00 203F 
00 F 

00 2560F 

00 800K 


1 

Decimal Order Inpi 


it 




o 


kl 202F 
LI 39F 


by 2 % from 3 






1 


L^ 202F 
kO 202F 




Compute sum check 




2 


F5 OL 
kO OL 








3 


LO 10L 
32 OL 








k 


00 IF 
L5 39F 


by 65 from 8 






5 


86 11F 










00 s4 


*>y 7 


Store Routine M-k on 




6 


F5 ^L 
42 ^L 




drum 




7 


F5 5L 
40 5L 








8 


LO 11L 
32 kL 








9 


26 999F 
00 F 




Read in another part 
of program 




10 


Nl 202F 
LI 202F 








n 


06 11F 

00 i6ksk 








12 


kl OF 
LI 39F ■ 


by 4 5 from 5 
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15 


LA OF 
kO OF 






Compute sum check 




Ik 


F5 12L 
k2 12L 










15 


LO 55L 
52 12L 










16 


■ 00 IF 












L5 59F 


by 18 5 


from 20 






17 


86 11F 












oo i6ksk 


by 195.' 




Store routine on drum 




18 


F5 16L 
ko i6l 










19 


F5 17L 
kO 1TL 










20 


LQ 5^L 
52 l6L 










21 


00 IF 
L5 OF 






Store sum check 




22 


86 11F 
00 225S4 






on drum 




25 


kl 59F 












LI 5F 


"by 25| 


from 26 






2k 


Lk 59F 
kO 59F 






Compute sum check 




25 


F5 25L 
k2 25L 










26 


LO 55L 
52 25L 










27 


00 IF 












L5 5F 


by 29^ 


from 51 






28 


86 11F 
00 226sk 


by 50 




Store routine on drum 




29 


F5 27L 
k2 2JL 
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30 


F5 28L 
40 28L 








31 


LO 36L 
32 27L 








32 


26 999F 
00 F 




Read in another part 
of program 




33 


Nl OF 
LI 100F 








3\ 


06 11F 
00 225SV 








35 


Nl 39F 
LI 39F 








36 


06 11F 
00 261S4 








37 


4l ll4F 
LI 82F 


by 395 fro 11 **° 






38 


L4 ll4F 

4o ii4f 


■ 






39 


F5 37L 
40 kjL 








ko 


LO kjL 
32 37L 








.'.hi 


00 IF 
L5 82F 


by ky 9 from 45 






k2 


86 11F 
00 26is4 


by 44 


Store routine on drum 




43 


F5 4ix 
40 4lL 








44 


F5 42L 
42 42L 








45 


LO 48L 
32 4lL 








46 


20 63F 
OOF 
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LOCATION 
^7 



3 
k 

5 
6 

7 
8 

9 
10 

ll 

12 



ORDER 

Nl ll^F 

Ll 114F 

06 IIP 

00 29^S4 

00 39K 

Library Routine M-4 

26 800N 



NOTES 



PAGE k 
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00 5K 
50 S3 
50 L 
26 77L 
.40 OP 
L3 4F 
52 3L 
FF 27F 
L3 OF 
32 5L 
L5 30L 
26 6L 
L5 31L 
kO 19L 
L5 3F 
LO 28L 
k2 29L 
00 20F 
k6 19L 
81 kF 
00 20F 
k6 18L 
kl 2F 
85 11F 
00 Sk 
32 12L 
kO (3^L) 



"by 60 5 from 73 



from 3 



from 4 
from 5 



Change OFF OFF in word 129 of 
M-4 to FF 28F FF 28F 



Read in data 



Test to see if enough elements 
have been read in 
Eigenvalues or 
Eigenvalues and Eigenvectors? 



Read in print 
parameter 



Read eigenvalue 
routine off drum. 
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ORDER 
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13 


iA- 2F 
^0 2F 






Ik 


F5 12L 
kO 12L 






15 


Fp XJLXi 
40 11L 






16 


LO 32L 
36 11L 






17 


L3 2F 




Has routine "been read 




32 18L 




off drum correctly? 


18 


FF (0)F 
92 I67F 


"by 10 




19 


50 ( )F 
50 19L 


by 6, 8 




20 


26 3te 

kl 2F 




Perform computations 


21 


85 11F 


from 26 


■ 




00 i64s4 


by 25 


Read next routine 


22 


32 22L 




off drum 




kO 3kh 


"by 2k 




23 


'lA 2F 
kO 2F 


. 




2^ 


F5 22L 
kO 22L 






25 


F5 21L 
kO 21L 






26 


LO 33L 








36 21L 






27 


L3 2F 
36 3^L 






28 


FF 30F 
00 IF 






29 


80 F 








00 (11 )F 


by 7 
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30 


JO F 
50 19L 






31 


50 F 
50 19L 






32 


05 IIF 
00 l64s4 






33 


05 HF 
00 226s4 






3* 


4l 5F 
L5 l&i 


"by 22| from 27 




35 


46 42L 
46 48l 




Print parameters. 


36 


L5 3BL 
42 4lL 






37 


L5 29L 
50 29L 






38 


74 29L 








S5 S3 




. 


39 


10 IF 
L4 38L 






40 


40 6F 








92 131F 


from 57 




in 


92 515F 








L5 ( )F 


by 55 




42 


50 ( )F 
«P0 42L 


by 35 




43 


26 77L 
LI 19L 


■ 


Enter print routine 


44 


36 54L 
L5 5F 




■ 


^5 


L4 6f 
42 4tl 






46 


92 131F 








41 7F 
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vr 


92 131F 


from 53 






L5 ( )F 


by 45,50 




48 


50 ( )F 

50 48L 


by 35 




49 


26 *nx, 
L5 ^7L 






50 


L4 29L 
42 kjL 






51 


92 515F 
F5 7F 






52 


kO JF 
LO 29L 






53 


36 47L 
92 131F 






54 


F5 4lL 
F4 5F 


from 44 




55 


42 4lL 
F5 5F 






56 


ko 5F 
LO 29L 






57 


32 40L 
92 551P 






58 


20 58L 




Stops prepare to "begin 




kl 2F 




new problem 


59 


85 11F 


from 64 






00 226S4 


by 63 




60 


32 60L 








40 OL 


by 62 




61 


L4 2F 
40 2F 






62 


F5 60L 
4o 60L 






63 


F5 59L 
40 59L 






6* 


LO 75L 
36 59L 
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ORDER 
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65 


L3 2F 
32 66l 






66 


FF 31F 








kl 3F 


from 65 




61 


85 11F 


from 72 






00 26lS^ 


by 71 




68 


36 68l 




■ ■ . 




40 77L 


by 70 




69 


Lk 3F 
kO 3F 






70 


F5 68l 






•' 


k2 6% 






71 


F5 67L 








ko 67L 






72 


■LO76L 




' .-' ... ■ - 




36 67L 






73 


L3 3F 

36 5F 






7^ 


FF 32F . 


■ 






00 F 


■ 




75 


05 11F 
00 261SV 




• 


76 


05 lli 1 




'.■'. 




00 29^34 




f 




00 82K 


, 


... 


Library Routine P~2 52 






26 812N 








00 82K 






Library Routine (Re\ 

1 


ised N-2) 




■ 00 900K 




Library Routine X-7 1^2 


Change FF order here to FF 26F 




26 837N 




. 


■ 


. 




■.■'■"' 
. ' " 








- - 
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